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An  anode  support  tubular  solid  oxide  fuel  cell  (SOFC)  is  fabricated  and  the  dependence  of  its  polarization 
resistance  on  anode  microstructural  parameters  is  investigated  by  means  of  stereology  and  concept  of  con¬ 
tiguity  (c-c)  theory.  Nickel  yttria-stabilized  zirconia  (Ni-YSZ)  anode  supported  cell  with  YSZ  electrolyte, 
lanthanum-strontium-manganite  (LSM)-YSZ  composite  cathode,  and  LSM  cathode  layers  is  fabricated 
by  dip  coating.  Submicrometer  resolution  images  of  anode  microstructure  are  successfully  obtained  by 
low  voltage  SEM-EDX  and  quantified  by  stereological  analysis.  Cell  voltage  measurements  and  impedance 
spectroscopy  are  performed  at  temperatures  of  650  and  750  °C  with  hydrogen  and  nitrogen  mixture  gas  as 
a  fuel.  A  quantitative  relationship  between  polarization  resistance  and  microstructural  parameters  such 
as  circularity,  three-phase  boundary  length,  contiguity,  etc.  is  investigated  using  the  concept  of  contigu¬ 
ity  (c-c)  theory.  The  effectiveness  of  correlating  polarization  resistance  of  anode  supported  tubular  SOFC 
using  stereology  and  c-c  theory  is  evaluated. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  oxide  fuel  cell  (SOFC)  is  one  of  the  most  promising  energy 
conversion  technologies  due  to  its  high  efficiency  and  fuel  flexibil¬ 
ity  [1].  In  order  to  reduce  manufacturing  cost  and  to  achieve  high 
reliability,  recent  activities  are  shifting  to  operate  SOFCs  at  interme¬ 
diate  temperatures.  Anode  supported  SOFCs  have  attracted  much 
research  interests  since  they  can  reduce  the  electrolyte  thickness 
as  well  as  the  electrolyte  resistance  at  low  temperature  operation 
[2-4].  Great  efforts  have  been  paid  for  developing  materials  with 
high  conductivities  and  electrodes  with  optimal  microstructures  to 
improve  the  electrode  polarization  characteristics  at  low  tempera¬ 
ture  operation. 

Porous  cermets  (such  as  Ni-YSZ  composite)  are  widely  used 
as  anode  materials  in  order  to  match  the  thermal  expansion 
coefficient  with  that  of  the  electrolyte  and  to  extend  effective 
reaction  sites  (three-phase  boundary,  TPB)  [5].  It  is  well  known 
that  anode  microstructure  has  significant  impact  on  cell  perfor¬ 
mance  and  durability  [6-10].  To  date,  great  efforts  have  been 
made  to  improve  anode  performance  by  modifying  the  anode 
microstructures.  For  example,  effects  of  morphology  improve¬ 
ment  [9,10],  dependence  on  starting  materials  and  fabrication 
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methods  [11-13]  have  been  reported.  Recently,  direct  measure¬ 
ments  of  SOFC  electrode  microstructures  have  been  reported 
using  focused  ion  beam  scanning  electron  microscopy  (FIB- 
SEM)  [14]  and  X-ray  computed  tomography  (XCT)  [15].  Through 
these  direct  measurements,  important  microstructural  parame¬ 
ters  such  as  three-phase  boundary  length  and  tortuosity  factor 
can  be  obtained.  However,  huge  efforts  are  required  for  obtain¬ 
ing  large  amount  of  data  volume  with  very  high  resolution  for 
statistical  processing.  A  methodology  that  can  quantify  microstruc¬ 
tural  parameters  from  two-dimensional  microscopy  images  is 
preferred  especially  from  the  designing  point  of  view.  The 
microstructural  parameters  such  as  phase  contiguity  and  effec¬ 
tive  conductivities  can  be  evaluated  from  two-dimensional  images 
using  stereological  theory  [16]  and  concept  of  contiguity  the¬ 
ory  [17,18].  Simwonis  et  al.  [19]  and  Lee  et  al.  [20-23]  have 
applied  stereological  theory  to  Ni-YSZ  cermets  and  investigated 
the  contiguity  of  the  phases.  Effects  of  Ni  agglomeration,  Ni 
contents  or  compaction  pressure  on  the  effective  conductiv¬ 
ity  of  the  Ni-YSZ  cermet  are  reported.  However,  such  approach 
for  the  electrodes  of  the  real  cell  configuration  is  still  limited. 
And  more  importantly,  the  quantitative  relationship  between  the 
microstructural  parameters  and  the  polarization  resistances  is 
not  fully  investigated.  Thus,  quantitative  investigation  between 
anode  microstructural  parameters  and  polarization  characteris¬ 
tics  for  the  real  anode  supported  cell  configurations  is  strongly 
demanded. 


0378-7753 /$  -  see  front  matter  ©  2009  Elsevier  B.V.  All  rights  reserved. 
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Nomenclature 

a  electrode  total  thickness  (|xm) 

Aj  average  area  of  phase  i  ( |xm2 ) 

Q  contiguity  of  phase  i 

CPE  constant  phase  element 

dj  equivalent  diameter  of  phase  i  ( pm) 

F  Faraday  constant:  9.649  x  104  (C  mol-1 ) 

i  current  density  (A  cm-2 ) 

li  mean  intercept  length  of  phase  i  (pm) 

Ltpb  three-phase  boundary  length  per  unit  volume 

(pm  pm-3) 

Nj_j  number  of  contact  points  of  phases  i  and  j  per  unit 
length  (pm-1) 

Ntpb  number  of  intersection  points  of  three  phases  per 
unit  area  (pm-2) 

Pi  average  perimeter  of  phase  i  (pm) 

R  resistance  (Q  m) 

R0  gas  constant:  8.314  (J  mol-1  K-1 ) 

Ri  circularity  of  phase  i 

Si  total  interfacial  area  of  phase  i  per  unit  volume 

(pm2  pm-3) 

S/_j  interfacial  area  of  phases  i  and  j  per  unit  volume 
(pm2  pm-3) 

T  temperature  (K) 

Vi  volume  fraction  of  phase  I 

Greek  symbols 
s  porosity  (=Vpore) 

r  Thiele  modulus 

a  conductivity  ( S  m- 1 ) 

Subscripts  and  superscripts 

eff  effective  value  in  multi  phase  media 

el  electron 

io  ion 

h  high  frequency 

1  low  frequency 

ohm  ohmic 

i  phase  i  (Ni,  YSZ-pore) 

Ni  nickel  phase 

Pore  pore  phase 

YSZ  YSZ  phase 

TPB  three-phase  boundary 


In  the  present  work,  microstructure  of  the  tubular  anode  sup¬ 
ported  SOFC  is  quantified  and  the  effects  of  anode  microstructure 
on  its  polarization  resistances  are  evaluated.  To  obtain  cells  with 
different  anode  microstructures,  two  different  NiO-YSZ  starting 
materials  are  used  and  the  amounts  of  poreformer  were  varied. 
Then,  electrolyte  and  cathode  layers  were  dip  coated  onto  the  anode 
supports.  I-V  curves  and  Nyquist  plots  were  obtained  at  tempera¬ 
tures  of  650  and  750  °C  using  hydrogen  and  nitrogen  mixture  gas 


as  a  fuel.  Reduced  anode  microstructures  were  imaged  using  low 
voltage  scanning  electron  microscopy  (SEM)  with  energy  disper¬ 
sive  X-ray  spectroscopy  (EDX).  Microstructural  parameters  such 
as  three-phase  boundary  (TPB)  length  and  contiguity  are  quanti¬ 
fied  by  means  of  stereology  [16].  The  effects  of  starting  material 
and  poreformers  on  microstructure  are  evaluated.  Finally,  the  rela¬ 
tionship  between  the  anode  microstructural  parameters  and  the 
polarizations  is  investigated  using  concept  of  contiguity  theory 
[17-19]. 

2.  Cell  fabrication 

Anode  supported  cells  used  in  the  present  study  consist  of  Ni- 
YSZ  anode  support  (thickness  500  jxm),  YSZ  electrolyte  (20  |xm), 
LSM-YSZ  composite  (10|xm)  and  LSM  cathode  current  collector 
(10  [Jim)  layers.  NiO-YSZ  (60/40  vol%)  powders  (AGC  Seimi  Chemical 
Co.,  Ltd.)  with  binder  (methylcellulose)  were  mixed  and  extruded 
into  a  tubular  shape  (Kankyo  Ceramics  Research  Co.  Ltd.,  Japan). 
Since  sintering  shrinkage  rate  is  sensitive  to  the  amount  of  pore- 
formers,  anode  pre-sintering  temperature  was  chosen  differently 
(1100  or  1200  °C)  for  the  cells  depending  on  the  anode  microstruc¬ 
tures  so  that  every  cell  has  nearly  the  same  electrolyte  thickness. 
Then,  YSZ  slurry  was  dip  coated  on  the  anode  support,  and  co¬ 
sintered  at  1400  °C  for  3h.  The  LSM-YSZ  (50/50  wt%)  composite 
cathode  layer  and  LSM  cathode  current  collector  layer  were  dip 
coated  and  sintered  at  1200  and  1100  °C  for  3  h,  respectively.  Finally, 
Ag  wire  was  connected  to  the  cell  by  Ag  paste.  Anode  side  cur¬ 
rent  collection  wire  was  connected  to  one  edge  of  the  anode 
support.  After  connecting  the  cell  with  the  fuel  supply  tube,  the 
sealant  (Dura  Seal  1529  UHT,  Cotronics  Corp.)  was  painted  over 
the  anode  with  the  current  collection  wire  for  perfect  sealing. 
Fig.  1  shows  the  schematic  of  the  fabricated  anode  supported 
cell. 

To  change  the  anode  porosity,  the  amount  of  acrylic  poreformer 
(1.5  [Jim,  Soken  Chemical  &  Engineering  Co.,  Ltd.)  was  varied  as 
0  wt%,  9  wt%  and  17  wt%.  For  the  starting  materials,  co-precipitated 
or  mechanically  mixed  NiO-YSZ  powders  were  employed.  The 
mean  grain  size  of  the  mechanically  mixed  powders  was  around 
1  [Jim.  Primary  particle  size  of  the  co-precipitated  powders  is 
30-50  nm  and  they  were  uniformly  distributed  inside  the  sec¬ 
ondary  grains  around  1  jxm.  Table  1  shows  the  specifications  of  the 
cells  tested  in  this  study.  Fig.  2  illustrates  the  cross-section  of  the 
fabricated  co-precipitated  cell  with  17wt%  poreformer  (C-17wt% 
cell). 

The  cells  were  connected  to  the  fuel  supply  tube  with  adhesion 
bond  (904  Zirconia,  Cotronics  Corp.)  and  sealant  (Dura  Seal  1529 
UHT,  Cotronics  Corp.).  Ag  wires  were  connected  to  the  electronic 
load  controller  (Kikusui  Electronics  Corp.,  PLZ664)  and  impedance 
meter  (Kikusui  Electronics  Corp.,  KFM2150). 

The  anode  was  reduced  in  H2  20  seem  and  N2  10  seem  at  750  °C 
for  10  h.  The  ramp  rate  was  10°Cmin-1  in  all  heating  and  cool¬ 
ing  processes.  Hydrogen/nitrogen  flow  rates  were  set  as  20/10  or 
2/28  seem  fixing  the  total  flow  rate  at  30  seem.  Air  flow  rate  was 
fixed  as  20  slm.  Then,  I-V  curves  and  Nyquist  plots  were  acquired 
at  temperatures  of  650  and  750  °C. 


(b) 

Sealant 

Fuel  ;  CfH  ■ifi/if/iiiil  t 

K 


Fig.  1.  Schematic  of  the  fabricated  cell. 
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Table  1 

Fabricated  anode  samples. 


Abbreviation 

Starting  NiO/YSZ  powder 

NiO/YSZ  ratio  (wt%) 

Amount  of  poreformer  (wt%) 

Pre-sintering  temperature  (°C) 

Co-sintering  temperature  (°C) 

C-0  wt% 

Co-precipitation 

60/40 

0 

1200 

1400 

C-9  wt% 

Co-precipitation 

60/40 

9 

1100 

1400 

C-17  wt% 

Co-precipitation 

60/40 

17 

1100 

1400 

M-0  wt% 

Mechanical  mixing 

60/40 

0 

1200 

1400 

M-9  wt% 

Mechanical  mixing 

60/40 

9 

1100 

1400 

M-17  wt% 

Mechanical  mixing 

60/40 

17 

1100 

1400 

Fig.  2.  Cross-section  of  the  fabricated  cell  (C-17  wt%). 


3.  Microstructure  quantification 

3  A.  Image  processing 

Samples  were  infiltrated  with  epoxy  resin  (Marumoto  struers 
K.K.)  under  vacuum.  Cured  samples  were  polished  by  Ar-ion 
beam  cross-section  polisher  (JEOL  Ltd.,  SM-09010).  This  technique 
provided  less  damage  and  smoother  cross-section  than  the  dia¬ 
mond  slurry  polishing.  Then,  microstructure  of  the  cross-section 
was  successfully  imaged  by  the  energy  selected  backscattered 
(ESB)  electron  detector  and  energy  dispersive  X-ray  spectroscopy 
(EDX;  Thermo  Electron,  NSS300)  in  low  voltage  scanning  elec¬ 
tron  microscopy  (SEM;  Carl  Zeiss,  ULTRA55).  The  size  of  the 
images  is  30.3  p,m  x  40.5  f±m  with  768  x  1024  pixels.  As  shown  in 
Fig.  3,  ESB  and  EDX  of  low  voltage  SEM  provided  images  with 
submicron  resolution.  Co-precipitated  cells  show  finer  and  more 
complicated  structures  than  the  mechanically  mixed  cells.  For 
image  processing,  a  weighted  differential  filter  was  applied  to 
sharpen  grain  boundaries.  Pore  was  first  distinguished  from  the 
ESB  images,  and  then  solid  Ni  and  YSZ  phases  were  distinguished 


using  the  EDX  images.  Fig.  4  shows  the  processed  images  of  the 
six  samples  listed  in  Table  1.  As  can  be  seen  from  the  figure, 
Ni  (white),  YSZ  (gray)  and  pore  (black)  phases  are  clearly  distin¬ 
guished. 

From  these  processed  images,  volumetric  fraction  Vit  area  Ait 
perimeter  Pit  mean  intercept  length  /,  of  each  phase,  number  of 
contact  points  of  i  and  j  phases  per  unit  length  Nj_j,  number  of 
intersection  points  of  three  phases  per  unit  area  NTPB  are  obtained. 
Diameter  and  circularity  Rt  of  each  phase  are  calculated  as: 


(1) 


Ri  = 


fL 

4 ttAj  ’ 


(2) 


Note  that  circularity  Rt  takes  larger  values  as  the  shape  becomes 
more  complicated. 


Fig.  3.  SEM-EDX  images,  red:  Ni,  green:  YSZ,  black:  pore,  (a)  C-0  wt%  cell  and  (b)  M-0  wt%  cell.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is 
referred  to  the  web  version  of  the  article.) 
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□  Ni  □  YSZ  ■  Pore  □  Ni  DYSZ  ■  Pore 


Fig.  4.  Processed  cross-sectional  images  of  the  anodes.  Co-precipitated  (a)-(c),  and  mechanically  mixed  (d)-(f)  anodes. 


3.2.  Stereology 

According  to  the  theory  of  stereology  [16],  the  microstruc¬ 
ture  information  of  two-dimensional  images  can  be  transformed 
to  three-dimensional  information.  Interfacial  area  S/_j,  which  is 
the  neck  area  between  i  and  j  phases  per  unit  volume,  and  TPB 
length  Ltpb  per  unit  volume  in  an  isotropic  structure  can  be  written 


as: 

SH  =  2  Nh, 

(3) 

Ltpb  =  2NTPB. 

(4) 

Contiguity  Q  [17-20]  is  defined  as  the  fraction  of  the  interfacial 
area  shared  with  the  same  phase  to  the  total  surface  area  of  the 
corresponding  phase.  Simiwonis  et  al.  [19]  extended  the  theoreti¬ 
cal  approach  for  the  equiaxed  two-phase  solid  to  the  Ni-YSZ-pore 


three-phase  problem  and  obtained  following  expression  for  conti¬ 
guity: 


Q  = 


Sj 

Si  +  Sj  +  sk 


(5) 


where  5Z  =  2Sj_j  +  St_j  +  St_k  is  the  total  interfacial  area  of  phase  i. 
Volume  to  surface  ratio  is  related  to  the  intercept  length  //  by: 


Then,  Eq.  (5)  can  be  rewritten  in  terms  of  volume  fraction  V/  and 
mean  intercept  length  as: 


Viljlk 

Viljk  +  Vjkli  +  VMj- 


(7) 


•v 

o 

n> 

i 

o' 


Fig.  5.  Performance  curves  for  H2  =  20  seem  at  750  °C.  (a)  Co-precipitated  and  (b)  mechanically  mixed  anodes. 
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Fig.  6.  Performance  curves  for  H2  =  2  seem  at  750  °C.  (a)  Co-precipitated  and  (b)  mechanically  mixed  anodes. 
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Fig.  7.  Impedance  spectra  at  i  =  0.16  A  cm-2,  (a)  C-0  wt%  and  (b)  C-9  wt%  cells. 


Thus,  contiguity  can  be  easily  obtained  from  the  two- 
dimensional  images  through  Eq.  (7).  The  effective  conductivity  will 
be  related  to  contiguity  Q  in  the  later  section. 

4.  Results  and  discussion 

4.1  Cell  testing 


shows  the  obtained  R h  and  R\  of  the  cells.  C-17  wt%  cell  and  M- 
17  wt%  cell  data  are  not  plotted  in  the  figure  because  of  the  leakage 
problem  described  above.  As  shown  in  Fig.  9,  co-precipitated  cells 
show  smaller  Rh  and  R\  than  the  mechanically  mixed  cells.  Addi¬ 
tion  of  poreformers  drastically  reduces  R\  but  increases  Rh.  Voltage 
drop  due  to  fuel  shortage  was  apparent  at  H2  =  2sccm,  T=650°C 
case  around  i  =  0.1  A  cm-2.  This  may  be  the  reason  for  the  large  R\ 
obtained  for  the  M-9  wt%  cell  at  H2  =  2  seem,  T  =  650  °C  in  Fig.  9(a). 


A  typical  example  of  the  obtained  voltage  and  power  density 
curves  at  750  °C  are  shown  in  Figs.  5  and  6.  As  shown  in  Fig.  5, 
co-precipitated  cells  shows  slightly  better  performance  than  the 
mechanically  mixed  cells  at  H2  =  20  seem,  and  the  maximum  power 
density  (MPD)  was  over  0.4  W cm-2  for  the  C-0wt%  cell.  C-17wt% 
and  M-17  wt%  cells  showed  low  open  circuit  voltage  (OCV).  From 
the  SEM  images,  it  is  observed  that  the  electrolyte  layers  of  C- 
17wt%  and  M-17wt%  cells  were  not  dense  enough,  which  lead  to 
imperfect  sealing.  Fig.  6  shows  the  results  for  the  H2  =  2  seem  case. 
Co-precipitated  cells  again  showed  better  performance  than  the 
mechanically  mixed  cells.  Especially,  C-9wt%  cell  exhibits  MPD  of 
0.15  W  cm-2  and  the  maximum  fuel  utilization  exceeded  85%.  Effect 
of  poreformers  was  different  depending  on  the  starting  materials 
and  the  hydrogen  flow  rates.  For  the  co-precipitated  cells,  C-0  wt% 
cell  performed  better  for  the  H2  =  20  seem  case,  while  C-9  wt%  cell 
was  better  at  H2  =  2  seem.  On  the  other  hand,  for  the  mechanically 
mixed  cells,  M-9  wt%  cell  performed  favorably  in  H2  =  20  seem  case, 
while  M-0  wt%  cell  was  slightly  better  at  H2  =  2  seem. 

To  separate  ohmic,  activation,  and  concentration  polarizations, 
impedance  spectroscopy  data  were  obtained  at  i  =  0.16  A  cm-2. 
Fig.  7  shows  a  typical  example  of  the  obtained  Nyquist  plot.  Two 
semi-circles  were  least-square  fitted  to  obtain  ohmic  (R0 hm),  high 
frequency  ( Rh )  and  low  frequency  {R\)  resistances.  The  equivalent 
electrical  circuit  is  shown  in  Fig.  8. 

The  sum  of  the  ohmic  resistances  of  the  electrolyte,  cathode  and 
anode  is  measured  as  R0 hm.  The  high  frequency  resistance  Rh  is  con¬ 
sidered  to  be  due  to  the  activation  overpotentials  of  the  anode  and 
the  cathode,  while  the  major  contribution  of  the  low  frequency 
resistance  R\  is  due  to  the  concentration  polarization  [24].  Fig.  9 


4.2.  Microstructure  quantification  results 

The  original  30.3  |xm  x  40.5  [xm  images  are  divided  into  four 
15.2  |xm  x  20.2  |xm  sub-regions  in  order  to  investigate  the  depen¬ 
dence  of  the  micro-structural  parameters  on  the  analyzed  area.  The 
parameters  directly  obtained  from  the  image  processing  are  listed 
in  Table  2.  The  mean  values  as  well  as  the  standard  deviations  of 
the  data  from  four  sub-regions  are  also  shown.  As  can  be  seen  from 
Table  2,  the  standard  deviations  of  the  four  images  are  quite  large. 
The  standard  deviations  of  the  parameters  related  to  the  pore  phase, 
e.g.  Apore  and  Ppore,  are  very  large,  while  those  for  the  solid  phases 
remain  relatively  small.  Table  3  shows  the  parameters  calculated 
from  Eqs.  (l)-(4)  and  (7)  using  the  data  listed  in  Table  2.  In  Fig.  10, 
equivalent  diameters  dx  are  plotted  against  porosity  s  (=Vp0re).  As 
the  data  show  quite  large  scattering,  data  from  four  sub-regions  are 
listed  in  the  figure  as  four  independent  points  with  a  same  sym¬ 
bol.  It  can  be  seen  from  the  figure  that  each  of  the  data  sets  for 
the  co-precipitated  cells  and  mechanically  mixed  cells  correlates 
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Fig.  8.  Equivalent  circuit  of  SOFC  anode. 


Table  2 

Microstructural  parameters  directly  obtained  from  image  processing. 


Cell 

Volume  ratio 

Area 

Perimeter 

Intercept  length 

Number  of  contact 
points  of  i  and  /  phases 

Interception  points 
of  three  phases 

Ntpb  (pan-2) 

VNi 

Wsz 

^/pore  (-e) 

ANi  (pan2) 

Aysz  (pan2; 

1  AporeCp^m2)  PNi  (p-m) 

PYSz  (pan) 

Ppore  (p-m) 

/Ni  (pan) 

Iysz  (p^m) 

Ipore  (p-m) 

I^Ni-pore 

(parr1) 

i^Ni-YSZ 

(parr1) 

NySZ-pore 

(pan-1) 

C-0  wt% 

Image  1 

0.375 

0.373 

0.252 

3.291 

1.143 

0.878 

8.343 

5.095 

4.990 

1.251 

0.757 

0.606 

0.413 

0.720 

1.235 

1.017 

Image  2 

0.351 

0.410 

0.239 

3.168 

1.654 

0.611 

8.973 

6.799 

3.557 

1.071 

0.820 

0.596 

0.426 

0.770 

1.235 

1.073 

Image  3 

0.350 

0.393 

0.257 

2.827 

1.746 

0.717 

7.926 

7.520 

4.063 

1.143 

0.742 

0.586 

0.441 

0.724 

1.261 

1.017 

Image  4 

0.322 

0.446 

0.232 

2.147 

2.442 

0.579 

5.826 

9.045 

3.311 

1.128 

0.912 

0.586 

0.346 

0.715 

1.233 

0.939 

Mean 

0.350 

0.405 

0.245 

2.858 

1.746 

0.696 

7.767 

7.115 

3.980 

1.148 

0.808 

0.593 

0.406 

0.732 

1.241 

1.011 

a 

0.014 

0.019 

0.009 

0.240 

0.325 

0.134 

0.527 

1.246 

0.727 

0.091 

0.042 

0.010 

0.014 

0.028 

0.015 

0.032 

<x(%) 

4.1% 

4.6% 

3.8% 

8.4% 

18.6% 

19.3% 

6.8% 

17.5% 

18.3% 

7.9% 

5.2% 

1.7% 

3.5% 

3.8% 

1.2% 

3.2% 

C-9  wt% 

Image  1 

0.292 

0.342 

0.366 

1.947 

1.039 

2.079 

5.854 

4.453 

6.890 

1.037 

0.769 

0.979 

0.358 

0.677 

1.049 

0.945 

Image  2 

0.277 

0.304 

0.419 

1.573 

0.881 

3.891 

4.772 

4.121 

11.578 

0.978 

0.727 

1.145 

0.376 

0.635 

1.064 

0.952 

Image  3 

0.292 

0.319 

0.389 

1.826 

1.304 

2.254 

5.142 

5.559 

6.559 

1.176 

0.794 

1.125 

0.349 

0.630 

0.966 

0.815 

Image  4 

0.307 

0.346 

0.346 

1.886 

1.167 

1.380 

5.281 

5.014 

4.862 

1.124 

0.760 

0.941 

0.360 

0.668 

1.098 

1.027 

Mean 

0.292 

0.328 

0.380 

1.808 

1.098 

2.401 

5.262 

4.787 

7.472 

1.079 

0.762 

1.048 

0.361 

0.653 

1.044 

0.935 

cr 

0.009 

0.019 

0.026 

0.191 

0.214 

1.000 

0.550 

0.753 

2.807 

0.102 

0.034 

0.091 

0.014 

0.026 

0.053 

0.077 

*(% ) 

3.0% 

5.8% 

6.9% 

10.6% 

19.5% 

41.6% 

10.5% 

15.7% 

37.6% 

9.4% 

4.4% 

8.7% 

3.8% 

4.0% 

5.1% 

8.3% 

C-17  wt% 

Image  1 

0.228 

0.220 

0.553 

1.589 

0.629 

14.124 

4.592 

2.956 

26.444 

1.131 

0.692 

1.662 

0.358 

0.427 

0.829 

0.740 

Image  2 

0.234 

0.256 

0.510 

1.563 

0.649 

6.015 

4.376 

2.787 

12.779 

1.098 

0.749 

1.515 

0.358 

0.446 

0.904 

0.769 

Image  3 

0.248 

0.287 

0.465 

1.858 

0.807 

4.193 

5.015 

3.314 

9.443 

1.212 

0.782 

1.426 

0.297 

0.503 

0.918 

0.740 

Image  4 

0.198 

0.279 

0.523 

1.192 

0.822 

4.864 

4.009 

3.715 

10.919 

0.897 

0.722 

1.433 

0.322 

0.489 

1.049 

0.756 

Mean 

0.227 

0.260 

0.513 

1.551 

0.727 

7.299 

4.498 

3.193 

14.896 

1.084 

0.736 

1.509 

0.334 

0.466 

0.925 

0.751 

cr 

0.010 

0.034 

0.044 

0.163 

0.098 

5.287 

0.325 

0.269 

9.009 

0.059 

0.046 

0.119 

0.035 

0.040 

0.048 

0.017 

or(%) 

4.6% 

13.0% 

8.6% 

10.5% 

13.4% 

72.4% 

7.2% 

8.4% 

60.5% 

5.4% 

6.2% 

7.9% 

10.5% 

8.5% 

5.2% 

2.3% 

M-0wt% 

Image  1 

0.408 

0.435 

0.157 

5.960 

1.992 

0.482 

13.246 

4.738 

3.047 

1.432 

1.337 

0.490 

0.521 

0.545 

0.696 

0.971 

Image  2 

0.396 

0.427 

0.177 

5.781 

2.427 

0.654 

12.280 

5.728 

3.857 

1.546 

1.416 

0.581 

0.527 

0.469 

0.736 

0.812 

Image  3 

0.419 

0.393 

0.188 

8.562 

1.676 

0.687 

15.876 

4.133 

3.829 

1.709 

1.321 

0.593 

0.498 

0.393 

0.759 

0.802 

Image  4 

0.315 

0.487 

0.198 

3.333 

2.871 

0.691 

8.153 

6.485 

3.932 

1.364 

1.476 

0.586 

0.479 

0.450 

0.870 

0.773 

Mean 

0.384 

0.436 

0.180 

5.909 

2.241 

0.628 

12.389 

5.271 

3.666 

1.513 

1.388 

0.563 

0.506 

0.464 

0.765 

0.839 

cr 

0.011 

0.022 

0.016 

1.556 

0.377 

0.110 

1.861 

0.805 

0.459 

0.139 

0.051 

0.056 

0.015 

0.076 

0.032 

0.095 

<x(%) 

3.0% 

5.1% 

8.7% 

26.3% 

16.8% 

17.6% 

15.0% 

15.3% 

12.5% 

9.2% 

3.6% 

10.0% 

3.0% 

16.3% 

4.2% 

11.3% 

M-9  wt% 

Image  1 

0.311 

0.383 

0.306 

3.976 

1.565 

1.880 

8.811 

4.717 

6.755 

1.426 

1.060 

0.894 

0.367 

0.453 

0.908 

0.675 

Image  2 

0.305 

0.328 

0.367 

4.455 

1.005 

2.562 

9.780 

3.317 

7.561 

1.433 

0.992 

1.174 

0.382 

0.426 

0.878 

0.760 

Image  3 

0.270 

0.355 

0.375 

3.314 

1.451 

2.397 

7.467 

4.310 

6.977 

1.293 

1.096 

1.109 

0.366 

0.355 

0.891 

0.649 

Image  4 

0.261 

0.349 

0.390 

2.670 

1.231 

3.416 

6.348 

3.769 

9.893 

1.286 

1.072 

1.141 

0.399 

0.358 

0.903 

0.668 

Mean 

0.287 

0.354 

0.360 

3.604 

1.313 

2.564 

8.102 

4.028 

7.797 

1.360 

1.055 

1.079 

0.379 

0.398 

0.895 

0.688 

cr 

0.022 

0.028 

0.038 

0.573 

0.296 

0.356 

1.161 

0.720 

0.416 

0.079 

0.053 

0.147 

0.009 

0.051 

0.015 

0.058 

<x(%) 

7.7% 

7.8% 

10.5% 

15.9% 

22.6% 

13.9% 

14.3% 

17.9% 

5.3% 

5.8% 

5.0% 

13.6% 

2.3% 

12.7% 

1.6% 

8.4% 

M-17  wt% 

Image  1 

0.286 

0.289 

0.425 

4.614 

0.963 

4.832 

9.404 

3.448 

13.098 

1.660 

0.899 

1.225 

0.384 

0.299 

0.924 

0.587 

Image  2 

0.292 

0.313 

0.396 

3.579 

1.020 

3.035 

8.562 

3.690 

9.300 

1.293 

0.920 

1.102 

0.446 

0.390 

0.959 

0.717 

Image  3 

0.302 

0.363 

0.335 

4.033 

1.236 

2.186 

9.203 

3.990 

7.370 

1.438 

1.079 

0.985 

0.365 

0.449 

0.937 

0.789 

Image  4 

0.302 

0.374 

0.324 

3.307 

1.383 

1.988 

7.380 

4.201 

7.242 

1.304 

1.108 

0.929 

0.420 

0.401 

0.954 

0.756 

Mean 

0.295 

0.335 

0.370 

3.883 

1.151 

3.010 

8.637 

3.832 

9.252 

1.424 

1.001 

1.060 

0.404 

0.385 

0.944 

0.712 

cr 

0.008 

0.038 

0.046 

0.519 

0.144 

1.351 

0.440 

0.272 

2.915 

0.185 

0.098 

0.120 

0.042 

0.076 

0.018 

0.102 

<r(%) 

2.9% 

11.2% 

12.4% 

13.4% 

12.5% 

44.9% 

5.1% 

7.1% 

31.5% 

13.0% 

9.8% 

11.3% 

10.4% 

19.6% 

1.9% 

14.4% 
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Fig.  9.  Rh  and  R\  from  impedance  spectroscopy  at  (a)  650  °C  and  (b)  750  °C. 
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Fig.  10.  Equivalent  diameter  d,-  against  porosity  s.  (a)  Ni,  (b)  YSZ  and  (c)  pore. 


Fig.  11.  Relationship  between  circularity  R,  and  equivalent  diameter  d,.  (a)  Ni,  (b) 
YSZ  and  (c)  pore. 
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Fig.  12.  Interfacial  area  between  different  phases  (a)  Ni-pore,  (b)  YSZ-Ni  and  (c) 
YSZ-pore. 


quite  well  with  porosity  even  porosity  shows  large  scattering.  This 
implies  that  the  variation  of  parameters  between  the  sub-regions  is 
mainly  due  to  the  variation  of  porosity.  Thus,  even  with  data  disper¬ 
sion,  it  is  expected  that  the  present  data  can  be  used  for  discussing 
the  differences  of  microstructural  parameters  between  the  cells. 
Fig.  11  shows  the  relationship  between  circularity  and  equivalent 
diameter  dz.  The  mechanically  mixed  cells  have  larger  Ni  diameter 
dNi  and  circularity  ftNi  than  the  co-precipitated  cells.  On  the  other 
hand,  co-precipitated  cells  have  larger  YSZ  circularity  ftysz  than 
the  mechanically  mixed  cells.  From  its  definition,  larger  circular¬ 
ity  corresponds  to  non-circular  shape  of  the  phase.  These  indicate 
that  YSZ  phase  of  the  co-precipitated  cells  is  more  complicated, 
while  Ni  phase  of  mechanically  mixed  cells  is  larger  and  somewhat 
agglomerated.  No  apparent  difference  between  co-precipitated  and 
mechanically  mixed  cells  is  observed  for  the  pore  circularity  Rpore  in 
Fig.  11(c).  Fig.  12  shows  the  interfacial  area  between  different  phases 
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Fig.  13.  TPB  length  LTPB  against  porosity  e. 

Sj_j.  Again,  each  of  the  co-precipitated  cell  and  mechanically  mixed 
cell  data  is  well  correlated  with  porosity.  Mechanically  mixed  cells 
show  slightly  larger  Ni-pore  interfacial  area,  while  co-precipitated 
cells  show  larger  YSZ-Ni  and  YSZ-pore  interfacial  areas.  From  above 
results,  it  can  be  concluded  that  good  connections  of  YSZ-YSZ  and 
YSZ-Ni  phases  are  achieved  in  the  co-precipitated  cells,  while  Ni 
phase  is  somewhat  isolated  and  agglomerated  within  the  pore  for 
the  mechanically  mixed  cells.  Fig.  13  shows  the  TPB  length  LTPB 
and  porosity  e  (=Vp0re).  Porosity  increases  but  TPB  length  decreases 
with  increasing  the  amount  of  poreformer.  Co-precipitated  cells 
have  longer  TPB  length  and  higher  porosity  than  the  mechanically 
mixed  cells.  The  TPB  length  ITPB  is  also  plotted  against  geometric 
average  of  the  equivalent  diameters  (dNi^Ysz^pore)1^3  in  Fig.  14.  It 
can  be  seen  that  all  LTPB  data  are  well  correlated  with  the  aver¬ 
age  equivalent  diameter  regardless  of  starting  material.  Thus,  the 
large  LTPB  of  co-precipitated  cells  can  be  attributed  mainly  to  the 
fine  microstructure,  especially  to  the  finer  Ni  phase  as  shown  in 
Fig.  10(a).  Even  the  cells  are  co-fired  with  the  YSZ  electrolyte  layer  at 
1400  °C,  it  is  considered  that  the  complicated  YSZ  network  and  large 
interfacial  area  between  YSZ  and  Ni  phases  of  the  co-precipitated 
cell  possibly  prevented  Ni  sintering. 

4.3.  Polarization  resistance  and  Microstructure 

Finally,  relationship  between  polarization  resistances  and 
microstructural  parameters  is  discussed.  According  to  Costamagna 
et  al.  [25,26],  Thiele  modulus  r  is  defined  as: 
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Fig.  14.  TPB  length  LTpB  against  geometric  average  of  equivalent  diameters 
(dNidYSzdpore)^3- 


Table  3 

Microstructural  parameters  calculated  from  Eqs.  (1 )— (4)  and  (7). 


Cell 

Diameter 

Circularity 

Interfacial  area 

TPB  length 

Contiguity 

dNi  (pm) 

dy sz  (pm) 

dpore  (pm) 

^Ni 

Rysz 

Ppore 

■^Ni-pore  (pm2) 

•^Ni-YSZ  (  Pm2 ) 

1  Sysz-pore  (pm2 

!)  Ltpb  (pm2  parr 

3)  cNi 

Cysz 

Cpore 

C-0  wt% 

Image  1 

2.047 

1.206 

1.058 

1.683 

1.808 
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Fig.  15.  Rh  plotted  against  R.H.S.  of  Eq.  (12).  (a)  650 °C  and  (b)  750 °C. 


where  z'o  is  the  lineal  exchange  current  density  per  unit  TPB  length, 
a  is  the  electrode  thickness,  and  o ®ff  are  the  effective  ionic 
and  electronic  conductivities,  respectively.  Thiele  modulus  repre¬ 
sents  the  ratio  between  the  characteristic  reactive  length  and  the 
actual  thickness  of  the  electrode  [26].  When  anode  thickness  is  large 
(T>  3),  effective  region  of  the  electrode  becomes  thin  and  will  be 
restricted  to  the  vicinity  of  electrode/electrolyte  interface.  Then,  the 
effective  electrode  conductance  becomes  inversely  proportional  to 
r  [25].  In  that  case,  anode  conductance  can  be  expressed  using  the 
linearized  Butler-Volmer  equation  as: 


1 

R 


iof'TPB^pdf  X 

KqI 


1 

T 


z'oH/tpb 


(9) 


From  the  concept  of  contiguity  (c-c)  theory  [17-19],  the  effective 
conductivity  of  a  multi-phase  material  (txeff)  can  be  expressed  as 
follows: 

aeff  =  'f2<TiCi\/i,  (10) 


where  Q  and  V\  are  conductivity,  contiguity  and  volume  frac¬ 
tion  of  phase  z,  respectively.  Since  ionic  conductivity  of  Ni  can  be 
neglected,  effective  ionic  conductivity  of  the  Ni-YSZ  cermet  can  be 
expressed  as: 

of0ff  =  ovszCyszVysz  •  ( H ) 

Substituting  Eq.  (11)  in  Eq.  (9)  and  neglecting  electronic  resistivity, 
the  anode  resistance  can  be  related  to  the  microstructural  param¬ 
eters  as  [27]: 


R 


PI°g)KoT  i 
ioFiTPB 


1 


\J  Cysz  Wszf  tpb 


(12) 


differences  of  anode  resistances  between  the  cells  can  be  attributed 
mainly  to  the  difference  of  the  microstructural  parameters  such  as 
contiguity,  volume  fraction  and  TPB  length.  Thus,  it  is  expected  that 
stereology  and  c-c  theory  is  effective  for  evaluating  the  quantita¬ 
tive  differences  of  microstructures  and  discussing  the  relationship 
between  microstructural  parameters  and  polarization  characteris¬ 
tics  of  anode  supported  SOFCs. 

5.  Conclusions 

In  the  present  study,  tubular  anode  supported  cells  with  Ni- 
YSZ  anode,  YSZ  electrolyte,  LSM-YSZ  composite  cathode,  and  LSM 
cathode  current  collector  layers  were  fabricated.  Amount  of  acrylic 
poreformer  and  starting  NiO-YSZ  powders  were  varied  to  obtain 
different  anode  microstructures.  Anode  cross-section  was  polished 
by  Ar-ion  beam  cross-section  polisher  and  successfully  imaged  by 
low  voltage  SEM-EDX.  Clear  images  of  low  voltage  SEM-EDX  made 
it  possible  to  distinguish  the  three  phases.  The  microstructural 
parameters  such  as  TPB  length,  contiguity,  etc.  are  quantified  using 
the  stereological  analysis.  Co-precipitated  cells  have  good  connec¬ 
tions  of  Ni-YSZ  and  YSZ-YSZ  phases,  which  lead  to  long  TPB  length 
with  high  porosity.  On  the  other  hand,  Ni  is  relatively  isolated  and 
agglomerated  within  the  pores  and  TPB  length  and  porosity  are 
small  in  the  mechanically  mixed  cells.  Impedance  spectroscopy  was 
performed  at  650  and  750  °C  with  hydrogen/nitrogen  mixed  gas 
as  a  fuel.  High  frequency  resistance  scaled  well  with  contiguity, 
volume  fraction  of  YSZ  and  TPB  length,  i.e.  1/a/CYszWsz^tpb-  It  is 
expected  that  stereology  and  c-c  theory  are  effective  for  evaluat¬ 
ing  the  microstructures  and  their  polarization  resistances  of  SOFC 
anodes. 
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